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ABSTRACT 

Since the discovery of quasi-periodic propagating oscillations with periods of order three to ten minutes in 
coronal loops with TRACE and SOHO/EYT (and later with STEREO/EUVl and HinodelEYS), they have been al- 
most universally interpreted as evidence for propagating slow-mode magnetoacoustic waves in the low plasma 
/3 coronal environment. Here we show that this interpretation is not unique, and that for coronal loops asso- 
ciated with plage regions (as opposed to sunspots), the presence of magneto-acoustic waves may not be the 
only cause for the observed quasi-periodicities. We focus instead on the ubiquitous, faint upflows at 50-150 
km/s that were recently discovered as blueward asymmetries of spectral line profiles in footpoint regions of 
coronal loops, and as faint disturbances propagating along coronal loops in EUV/X-ray imaging timeseries. 
These faint upflows are most likely driven from below, and have been associated with chromospheric jets that 
are (partially) rapidly heated to coronal temperatures at low heights. These two scenarios (waves vs. flows) 
are difficult to differentiate using only imaging data, but careful analysis of spectral line profiles indicates that 
faint upflows are likely responsible for some of the observed quasi-periodic oscillatory signals in the corona. 
We show that recent EIS measurements of intensity and velocity oscillations of coronal lines (which had pre- 
viously been interpreted as direct evidence for propagating waves) are actually accompanied by significant 
oscillations in the line width that are driven by a quasi-periodically varying component of emission in the 
blue wing of the line. This faint additional component of blue-shifted emission quasi-periodically modulates 
the peak intensity and line-centroid of a single Gaussian fit to the spectral profile with the same small ampli- 
tudes (respectively a few percent of background intensity, and a few km/s) that were previously used to infer 
the presence of slow mode magneto-acoustic waves. Our results indicate that it is possible that a significant 
fraction of the quasi-periodicities observed with coronal imagers and spectrographs that have previously been 
interpreted as propagating magnetoacoustic waves, are instead caused by these upflows. The different physical 
cause for coronal oscillations would significantly impact the prospects of successful coronal seismology using 
propagating disturbances in coronal loops. 

Subject headings: Sun: chromosphere - Sun: corona - Sun: oscillations - Sun: magnetic fields 



1. INTRODUCTION 

Using observations from instrum ents on the Sol ar and He- 
liospheric Observatory (SOHO; Fleck et al.|1995|l, the Tran- 



sition Region and Coronal Explorer (TRACE; |Handy et al. 



1999 ), twin STEREO SECCHI/EUVI imagers ( [Howard et aL 
2008), and the Extreme -ultraviolet Imaging S pectrograph 



(EIS; |Culhane et al.|2007| ) on Hinode ( |Kosugi et al.[2007| > the 
community has invested a great deal of effort in the identi- 
fication and analysis of low-amplitude wave-like phenomena 
seen in EUV coronal imaging and spectroscopy (see, e.g.,|de 



Moort el etalTiO OO, 2002b a; King et al.|2003i|Andries et aL 
2005 l |McIntosh et al.|2008] 1Ofman & W ang|2008[ Wang et al. 
2DU8l[ Mariska et al.||2008[ |Van Doorsselaere et al.||20 08a b; 
Marsh et al. 2009 ; Marsh & Walsh 2009 ; Wang et al. 2009b a ; 
Mariska & Muglach 2010 ] for a few of the most cited, and re- 



cent examples). Nakariakov & Verwichte (2005) provides an 
excellent overview of the techniques used and highlights the 
community's interest in isolating and characterizing coronal 
waves, in order to remotely sense the physical attributes of 
the outer solar atmosphere by studying propagation speeds, 
amplitudes, and phase relationships of the observed phenom- 
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ena. 

The broad community effort to probe the properties of 
the coronal plasma and magnetic field by attributing MHD 
wave properties to (quasi-) periodic propagating disturbances 
will be accelerated by the availability of considerably higher 
signal-to-noise, high cadence, high spatial resolution, multi- 
wavelength imaging provided by the Atmospheric Imaging 
Array on the Solar Dynamics Observatory (SDO). It is the 
availability of this new, complex, and rich data in concert with 
recent investigations of chromospheric-coronal coup ling ( |De| 
|Pontieu et aL|2009||McIntosh & De Pontieu|200"9b] l that mo- 
tivate the work presented and the plainly stated caveat emptor 
warning that it carries: "not everything periodic in the outer 
atmosphere is evidence of a wave". 

In some sense we are revisiting an issue that was trig- 
gered by the first observations of propagating disturbances 
with TRACE and SOHO/EYT, with early reports suggesting 
that fl ows may explain so me of the observed disturbances 
(e.g. Schrijver et al.|1999] l. However, the presence of coher- 
ent 3 minute oscillations emanating from sunspots (clearly 
associated with umbral oscillations), and the lack of spec- 
troscopic evidence for the required high-velocity (^50-100 
km/s) line-of-sight flows, quickly led to the dominant inter- 
pretat ion of propagating slow mode magneto-acoustic waves 
(e.g., |de Moortel et aLl |2000). This interpretation ignores 
the fact that many of the observed propagating disturbances 
(especially in coronal loops emanating from quiet Sun net- 
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Figure 1. TRACE 195A intensity (panel A), inferred single Gaussian fit parameters to the EIS spectroheliogram showing the peak intensity (panel B), (relative) 
Doppler velocity (panel C), Gaussian width (panel D) and the results of the 75-125km/s R-B analysis in Fe XIII 202A. On each of the panels we also show the 
pointing of the timeseries (vertical solid line), and a dashed box showing the maxima of the EIS pointing drift for the timeseries observation. Panels A and E of 
this figure are supported by movies in the online edition of the journal showing the pointing variation of EIS on the TRACE images and the complete range of the 
R-B analysis. 



work and active region plage) do not show evidence of sig- 
nificant quasi-periodic signals. More importantly, inspired by 
recent observations of high speed, but faint upflows in spec- 
tra o f the TR and co ronal footpoints of active region loops 
(e.g., Hara et al7||2008 1, and of pervasive propagating distur - 
bances in imaging data of coronal loops ( |Sakao et al.|2007) , 
we discovered that faint upflows are ubiquitous in footpoints 
of coronal loops in activ e regions (|De Pontieu et al.||2009] l 
and quiet Sun alike (Mcintos h & De Pontieu||200'9lif . Us- 
ing ///not/e/SOT-EIS-XRT data, we have suggested a direct 
link between faint blueward asymmetries in TR/coronal spec- 
tra at loop footpoints, propagating disturbances in images of 
coronal loops, a nd often quasi-periodically recurring chromo- 
spheric spicules ( |De Pontieu et al.|2009[|McIn tosh & De Pon- 
|tieu|2009b| . Armed with these new findings, we use EIS spec- 
tra (described in § 2) and re-analyze them with two novel anal- 
ysis techniques: asymmetry analysis (§ 3), and guided dou- 
ble fits using a genetic algorithm (§ 4). We use Monte Carlo 
simulations to show (§ 5) how the observed variations of line 
intensity, centroid and line widtrrl are fully compatible with 
the presence of a faint, strongly blueshifted component that 
quasi-periodically recurs. We discuss the interpretation and 
impact of our findings in § 6. We show that the complexities 
of the atmospheric coupling, while appearing to detract from 
"coronal-seismology", offers the potential for a deeper under- 
standing of the mass and energy flow between chromosphere 

4 Hereafter we will simple refer to the line width as the |erm fpr the Gaus- 
sian width (a) of the line profile - I(X, t) = 7 e" (A ~ Ao) /< 2<T o). 



and corona. 



2. OBSERVATIONS & ANALYSIS 



The primary data set investigated is a combination 
of TRACE and Hinode/ElS observations of active region 
AR10940 on February 1, 2007 . Some of this data was dis- 



cussed by Wang et al. (2009b|l. While the TRACE observa- 



tions, taken in the 195 A passband (at a cadence of 74s and 
exposure times of 65s) are impacted by cosmic ray hits they 
provide very useful context for the interpretation of the EIS 
observations. For both instruments we perform the usual post- 
processing: correcting for cosmic ray hits, hot pixels, detector 
bias, and dark current, and, for EIS, converting the data num- 
bers to physical intensities (in erg cm -2 s" 1 sr -1 A -1 ). In 
addition to the regular preparation we use the solarsoft routine 
eis_jitter.pro to compute the pointing jitter of EIS using 
the spacecraft housekeeping data. As a result we can reduce 
the impact of the pointing jitter of the spectroheliogram and 
timeseries observations, and co-align TRACE and EIS. 

There are two different types of EIS observations: the first 
is an 256" x 256" (West to East) spectroheliogram of the re- 
gion starting at 00:12UT with a 60s exposure time at each 
rasterstep. The second part is a "sit-and-stare" (within the 
limits of the pointing jitter) observation starting at 01:32UT 
that involves 350 steps of the 1" x 512" slit with 60s ex- 
posures and solar rotation tracking at a median pointing of 
-140" (±24), -5" (±256). Panels A through D of Fig. [^pro- 
vide the imaging context for the timeseries observation: Panel 
A shows the TRACE 195 A image at the start of the spectrohe- 
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Figure 2. Single Gaussian fit parameters to the EIS timeseries. From top to 
bottom: the peak intensity, (relative) Doppler velocity, width, and the results 
of the 75-125km/s R-B analysis in Fe XIII 202A. On each of the panels we 
also show the region (y = 60-65") along the slit chosen for further analysis. 

liogram observation with the EIS timeseries slit position (ver- 
tical thin line) and maxima of the pointing jitter drift shown as 
a dotted box; Panels B, C, and D show the peak line intensity, 
(relative) Doppler velocity and width that result from a single 
Gaussian fit to the spectra of the Fe XIII 202A lin e which is 
formed (und er equilibrium conditions) at ~4.2MK (Mazzotta 
et al.||1998) . In this case we choose to analyze the Fe XIII 
202 A line instead of the Fe XII 195 A line used by Wang et al. 
( 2009b j) because the latter is affected by spectroscopic blends 
(e.g., |Brown et al.]|2008| l. This means that single Gaussian 
spectral fitting and studies of profile as ymmetry are more re- 
liable for the Fe XIII 202A line (e .g.. [Mcintosh & De Pon- 
fieu|20 09b; Mci ntosh et al.|2010a| > than for the Fe XII 195 A 
line. The similar formation temperatures of both lines mean 
that we can r eliably compare our results with those of | Wang 
et al.| ( p009b[ ). We will also exploit the spectroheliogram and 
time series observatio ns in the Fe XIV 264 and 274A lines 
(formed at — 2.5MK, |Mazzotta et al.||1998| ) to demonstrate, 
using more spectrally clean lines with relatively high signal- 
to-noise (S/N), that the observed phenomena are not isolated 
to a narrow temperature range. Unfortunately, while we ob- 
serve similar diagnostics and variability in other clean, cooler 
spectral lines in these data (e.g., Si VII 275 A) they have lower 
S/N and do not make for the most appropriate demonstration 
and so their presentation is omitted here. 

2.1. Line Profile Asymmetry Analysis 

Following the descript i on of|De Pontieu et al.| (|2009[l and 
Mcintosh & De Pontieu (2009b) we perform a Red Minus 
Blue' (R-B) profile asymmetry analysis on spectral lines in 
the EIS data that are not significantly impacted by spectral 
blends in the relatively narrow (24 pixel) spectral windows 
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Figure 3. Five-minute running differences for the EIS timeseries shown in 
Fig. [2] From top to bottom: the peak intensity, (relative) Doppler velocity, 
widfh^and the results of the 75-125km/s R-B analysis in Fe XIII 202A. 

used, e.g., the Fe XIII 202A line. The R-B analysis involves 
several steps. First we fit a single Gaussian to the emission 
line profile to establish the line centroid. Once determined, 
we sum the amount of emission in narrow (^24 km s _1 wide) 
spectral regions symmetrically placed about the determined 
centroid. When summing we use a line profile that is inter- 
polated to a spectral resolution that is ten times finer than that 
of EIS. We then subtract the red and blue wing contributions 
of the interpolated profile to make a filtergram that samples a 
particular velocity range. A positive value of R-B indicates an 
asymmetry in the red wing of the line, which we can interpret 
as the signature of excess down-flowing material at a velocity 
of that order while, conversely, a negative value of R-B would 
indicate an excess of up-flowing material. Panel E of Fig 1 
shows the R-B diagnostic for Fe XIII 202 A at 75-125km/s and 
we are immediately drawn to the several locations where sig- 
nificant high velocity upflow signatures exist. Some of these 
underly the slit position of the time-series. Visual inspection 
of the TRACE movie accompanying panel A of Fig.[T]shows 
that this is a location of propagating intensity dist urbances (or 
so-called "blobs"), e xac tly like those discus sed in Mcintosh & 
|De Pontieu! (|2009b|) (or|Sakao et al.|2007] l. Indeed, the anal- 
ysis of I Wang et aLj (j2009b ) shows that the amplitude of the 
apparent motion of these blobs (deduced from motion along 
the loop in the plane of the sky) is of the same order as the 
velocity range that we have chosen for the R-B map shown. 



3. TIME-SERIES ANALYSIS 

To facilitate comparison, we select th e same portion of th e 
de-jittered timeseries as that chosen by Wang et al. ( 2009b ), 
i.e., from 03:06-06:00UT. Figure|2]shows the single Gaussian 
fit parameters for the Fe XIII 202A timeseries observations, 
from top to bottom showing the peak line intensity, (relative) 
Doppler velocity inferred from the shift of the profile centroid, 
the width and the 75-125km/s R-B analysis. We note the pres- 
ence of a large region of blue-wing asymmetry in this line at 
positions ~ y=-90" and around 60<y<110" that are consis- 
tent with locations identified in panel E of Fig. [T] The latter 
location is at the footpoint region of a loop fan that is clearly 
visible in the TRACE 195A image in Fig. [I] and on which 



4 



De Pontieu & Mcintosh 



FE XIII 202.M0A 




FE XII195.120A 



£ 5 

I 
a. 

i 

£ 



FE XIV 274.200A 



g 5 
1 




1-Feb-2O07 houfs artof midnight 

Figure 4. Extracted EIS timeseries for the Fe XIII 202A (top), Fe XII195A 
(middle) and Fe XIV 274A (bottom) emission lines from 03:06-06:00UT, for 
locations y = 60 — 65". In each panel we show the detrended percentage 
change in signal amplitude or velocity change in km/s. The green trace shows 
the temporal variance in the peak line intensity while the black, red, and 
blue traces show the changes in total profile intensity, Doppler velocity (line 
center shift), and line width respectively. Compare and contrast with the 
black dashed and dotted lines which show the 74-157km/s R-B analysis and 
detrended form (as for V and V n t). The zero lines for these are drawn at -2 
and -4 respectively and to stress our convention that values above that zero 
line represent an excess in the blue wing of the line. For clarity, the timeseries 
for intensity is offset by +4, velocity by +2 and the bottom linewidth by -2 
from the zero-point on the y-axis. 



blobs are seen to propagate in the TRACE and EIS data. The 
pointing jitter in the direction perpendicular to the slit most 
likely leads to the gradual, but significant decrease of the RB 
signal towards the end of the timeseries. We will discuss the 
effects of this pointing jitter further in § 6. 

A closer inspection of the peak intensity, the first and sec- 
ond moments of the Gaussian fit (i.e, Doppler velocity and 
line width) reveal episodes of small amplitude oscillations of 
a few percent in intensity and a few km/s amplitude for both 
the Doppler velocity and width. The oscillations of the peak 
intensity and Doppl er velocity were reported previously by 
Wang et al.| ( |2009b| for the blended Fe XII 195 A (cf. their 
Fig. 2 where, using running differences and detrending, they 
are highlighted beautifully). Here we highlight the (previ- 
ously unreported) existence of oscillations in the line width 
that are in phase with those of the peak intensity and veloc- 
ity, and that are accompanied by oscillatory changes in the R- 
B signal, with blueward asymmetries occurring in phase with 
line width increases. The oscillations in the line widths and R- 
B signal are weak, but visible in the 5-minute running differ- 
ences shown in Fig. [3] although with substantially lower con- 
trast than the lower order Doppler velocity and line intensity 



measurements^] The sign-flipping diagonal striations in the 
space-time plots, especially those around 60" < y < 150" of 
all four measurables are clearly related. As we will illustrate 
below, these line width and blueward asymmetry oscillations 
are very valuable diagnostics that provide strong physical con- 
straints on what is driving the low amplitude, quasi-periodic, 
signal in the data. 

The correlation between intensity, velocity, line width and 
blueward asymmetry is more clearly illustrated in Fig. [4] This 
figure shows an example oscillatory region from the time- 
series at the y =60-65" location as shown by the horizontal 
lines in Fig. [2j We have summed o ver the same region that 
was studied by [Wang et al. ( 2009b i in t heir Fig. A2. In the 
left panel of Fig. H we use the recip^Jof Wang et al. (2009b) 
and show the resulting, detrended variations in the peak in- 
tensity (black solid line), total intensity (here defined as the 
peak intensity times the line width; green solid line), Doppler 
velocity (red solid line), Gaussian line width (blue solid line), 
the R-B measure for velocities of 75-125km/s (black dashed 
line) and its detrended form (black dotted line). We repeat 
this process for the Fe XII 195 A line (middle) and the Fe XIV 
274A line (bottom), to illustrate that these variations are not 
isolated to the formation temperature of Fe XIII. Again, al- 
though not shown here, analysis of Si X 261 A and Si VII 
275A reveal a similar behavior, but are omitted here because 
those lines have significantly lower S/N. We also note that 
the presence of a slight blend (and/or slight gradient in the 
background emission) in the red wing of Fe XIII leads to oc- 
casional more slowly evolving excursions of the absolute RB 
measure towards the red. Detailed analysis suggests that the 
blend is likely caused by a much hotter line. The decent cor- 
relation between the RB asymmetry and the other moments of 
the Fe XIII, and with the RB measure in other lines (Fe XIV 
and Si X) indicates that the detrended RB measure is still a 
good measure of blueward excursions, especially since the 
redward excursions happen on longer timescales than the os- 
cillations we focus on here. 

Detailed analysis of the Fe XIII and Fe XIV panels in 
Fig. [4] show that the correlation between intensity, velocity, 
line width and blueward asymmetry is by no means perfect, 
but significant throughout the 3 hour long timeseries. Most 
of the line width peaks are associated with stronger blueward 
asymmetries and with blueward excursions of the line cen- 
troid, and intensity peaks. These correlations will be explored 
further in a statistical sense in § 4. We also note that the results 
of Fig. [4] are by no means unique: we have found many other 
locations in the same dataset that show such correlations. Two 
examples are shown in the appendix (Figs.[T][2]i. 

Detailed comparison of the time variations of these lines 
also shows that the oscillations of the intensity and Doppler 
velocity of Fe XIII 202 A and Fe XII 195 A have similar ampli- 
tudes and are very well correlated. While the correlation with 
Fe XIV 274A is not as clean, it is significant. This suggests 
that the driving mechanism for these oscillations acts over a 
wide range of temperatures. This does not seem compatible 
with a scenario in which sound waves propagate thro ugh an 
isothermal medium, as suggested for other datasets ( |Marsh &] 



5 The reduced contrast for higher order moments is understandable given 
the impact of photon noise on these measurements. This is illustrated in detail 
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Figure 5. Wavelet power spectra for the Fe XIII 2Q2A timeseries shown in the top panel of Fig 1 3 1 The solid white contours denote regions of the wavelet power 
spectrum of 95% significance and the cross-hatched region encloses the cone-ot-influence tor thepower spectrum. 
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Figure 6. Simulated timeseries of spectral line parameters from a forward 
model that simulates the effect of a faint and periodic secondary Gaussian 
component on the total intensity (full black line), peak intensity (green), line 
centroid (red) and line width (blue) from a single Gaussian fit to the total 
spectral line profile (which is the sum of two Gaussians reduced to EIS spec- 
tral resolution at 202A). The relative amplitude of the secondary component 
is 5% of the stationary background component, with a\ = a% = 30 km/s, 
and the second component is offset 8v = 50 km/s to the blue of the dominant 
component. We use the same convention as in Fig.p] in each panel we show 
the detrended percentage change in signal amplitude or velocity change in 
km/s. The black dashed and dotted lines show the 55-155km/s R-B analysis 
and detrended form (same recipe as for line centroid and width). The zero 
lines for these are drawn at -2 and -4 respectively and to stress our conven- 
tion that values above that zero line represent an excess in the blue wing of 
the line. For clarity, the timeseries for intensity is offset by +4, velocity by +2 
and the bottom linewidth by -2 from the zero-point on the y-axis. The input 
periodic signal has a period of P/2 = 6 minutes. The upper panel assumes 
that there is no photon noise, the middle and bottom panels have the same 
input periodic secondary component, but with photon noise added. The mid- 
dle panel is for S/N = 20, which is similar to the quality of a measurement 
of an individual EIS pixel in the dataset studied (l-Feb-2007), whereas the 
bottom panel is for S/N = 45, which is representative of the average over 
5 pixels that is shown in Fig. [4] The simulated profiles (with noise added) 
agree very well with the data ofFig.|4] both with respect to the average am- 
plitude of the variations, and the less-than-ideal correlations between some 
of the parameters. The agreement is remarkable when taking into account the 
fact that the Sun most likely does not produce a secondary component that is 
perfectly periodic and has a constant amplitude. 



Wals h|2009) l. 

The temporal behavior of these single Gaussian fit param- 
eters strongly suggests quasi-periodic behavior. The wavelet 



analysis (e.g., Torrence & Compo 1998) presented in Fig. [5] 
show, from left to right, the normalized wavelet power spectra 
of the intensity, doppler velocity, line width, and R-B asym- 
metry timeseries of Fig. [4] Outside the Cone-of-Infiuence (the 
white cross-hatched regions) and inside the 95% signal sig- 
nificance levels (solid white contours that are estimated us- 
ing a red noise background) we see that periodicities occur 
throughout the 3 hour long timeseries. These timeseries are 
best described as quasi-periodic with periods in a range from 
five to twelve mi nutes dominating - similar to those reported 
for Fe XII 195A ( |Wang et al.|2009bj) . This confirms the sim- 
ilarity in the oscillatory behavior of the Fe XIII and Fe XII 
lines. The quasi-periodicities in the R-B asymmetry suggests 
that the blueward asymmetries and associated upflows recur 
quasi-periodically. 

4. FORWARD MODELING 

What causes these oscillations? The close correlation be- 
tween blueward asymmetries and increased line width (as 
well as blueshifted line centroid and increased peak inten- 
sity) shown in Fig. [4] provides a strong clue. We use forward 
modeling to show that these observations are fully compati- 
ble with a scenario in which the quasi-periodic occurrence of 
a strongly blueshifted, but faint emission component causes a 
blueward asymmetry (R-B signal) at high velocity, and var- 
ious changes in the line parameters deduced from a single 
Gaussian fit to the spectral profile (as is common practice): 
slight increases in the line width, the position of the line cen- 
ter (the Doppler velocity), and the peak intensity of the line. 

To investigate this scenario, we use Monte Carlo simula- 
tions in which we construct a timeseries of simulated spectra 
with properties similar to the EIS observations. First, we cal- 
culate for each timestep the emergent spectrum assuming the 
spectral line emission is dominated by two different sources, 
both with Gaussian profiles (as a function of wavelength). The 
spectral line profile is given by the sum of two Gaussians: 

/(A, t) = I e-^- x ^^ 2 ^+aI cos 2 {2vt/P) e -^f/^l) 

(1) 

where Io, Ao, and ero are, respectively, the intensity, con- 
stant central wavelength, and line width of the steady back- 
ground component, and where a (<< 1) is the relative ampli- 
tude of the fainter second component compared to the back- 
ground component, Ai and <j\ are the center position and line 
width of the second component. The intensity of the sec- 
ond component is allowed to change with time (with a period 
P/2). This represents the recurring upflowing component that 
modifies the signal with time. 

To obtain ElS-like spectra, we convolve I(X,t) with a 
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Figure 7. Results of a parameter survey of the forward model in which we investigate the impact (on single Gaussian fit parameters) from a periodically recurring 
secondary component superimposed on a stationary background component. The properties of the two Gaussian components to the spectral line profile are varied 
over a wide range: 0.03 < a < 0.09 and Sv = 10 — 80 km/s, with fixed <tq = 30 km/s and u\ = 20 km/s. The secondary component is periodically recurring 
with a period P/2 = 6 minutes. The top row shows, as a function of a (amplitude) and Sv (velocity offset), the standard deviation over a one hour time range 
for the following fit parameters from a single Gaussian fit to the spectral line profile that is reduced to EIS spectral resolution: total intensity (defined here as 
peak intensity times line width), peak intensity, line centroid and line width. The bottom row shows the ratio between peak and total intensity, line centroid to 
line width, peak intensity to line centroid, and blueward asymmetry to peak intensity. The contours show the range of parameters that are "allowed" based on 
measurements with EIS (see text for details), with black and white contours respectively for upper and lower values of the allowable range. The hashed out 
regions show the non-allowed values. The bottom right panel shows the combination of all observational constraints, with only a small region for amplitudes 
of 4 — 8% and offset velocities of order 30-50 km/s. An accompanying movie (available in the online edition of the journal) illustrates that while the shape of 
the region of allowed parameters changes somewhat for different values of cr, the range remains locked around amplitudes of order 3-8% and velocities of 30-50 
km/s. 



Gaussian that has a (Gaussian) line width equal to the in- 
strumental broadening of EIS (22.93 mA). Next, we rebin 
the high resolution spectra to EIS resolution (with one EIS 
pixel equ al to 22 .93 mA). For each timestep we add photon 
noise t)) to each spectral pixel using random num- 

bers from a Gaussian distribution given by the IDL function 
randomn.pro. Next, we use gaussf it. pro to fit a single 
Gaussian to the emergent profile, and detrend the synthetic 
time series in a fashion identical to that used to prepare Fig. [4] 
Finally, we perform an R-B asymmetry analysis identical to 
that performed earlier on the line profiles observed with EIS 
(see § 2 for details). 

We repeat this recipe for a wide range of values for a, cr , 
<7i, and Sv[= c(\q — Ai)/Ao]; where c the speed of light]. An 
example timeseries for the peak intensity, total intensity, line 
centroid, line width and R-B asymmetry analysis for Sv = 50 
km/s, a — 0.05 and 0-0 = 0-1= 30 km/s is shown in Fig. [6] 
The top row shows, for this case, and when no noise is present, 
that the oscillations in peak and total intensity are of order 
a few percent, whereas the oscillations in line centroid and 
width are of order 1 km/s. All of these oscillations are in 
phase with each other, with the driver (the periodically recur- 



ring faint second component), and with oscillations in the R-B 
asymmetry: peaks in intensity occur at the same time as blue- 
ward excursions of the line centroid, the peaks in line width 
and the blueward asymmetries. To allow direct comparison 
we have produced plots that are identical in range in x and y 
to Fig. |4] Comparing the top row of Fig. [6] and Fig. |4] shows 
that generally the appearance of the "no noise" case is similar 
to what we observe with EIS. 

This correspondence (in a statistical sense) is even more 
striking when we compare Fig. [4] with the two other rows 
of Fig. [6] We have calculated this particular realization for 
three different levels of signal to noise (S/N): infinite (no 
noise), S/N=20, and S/N=45. Using the errors calculated by 
eis_prep.pro, we estimate that single pixels in the region of 
the observed oscillations (in the vicinity of y = 60 — 65") 
have a S/N of order 20. However, the oscillations shown in 
Fig.|4]were based on summing over 6 spatial pixels and thus 
have a larger S/N of order 45. The summing in this region 
was performed to al low direct comparison with the results of 
Wang et al. ( 2009b ) who suggested that the oscillations may 
be coherent in this region. We see that the S/N = 45 case 
shows an amazing correspondence (in a statistical sense) with 
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the oscillations observed with EIS in the Fe XIII 202 A line. 
This agreement is even more remarkable given the fact that 
the Sun most likely does not produce a perfectly periodic, 
constant amplitude secondary component. We also note that 
the increased noise-level has the strongest impact (not sur- 
prisingly) on the higher order moments of the profile: line 
centroid and width. It also has a significant effect on the R-B 
measure, and a somewhat smaller effect on the peak intensity. 
The increased noise leads to a correlation that is significantly 
worse than the "no-noise" case: random phase shifts between 
parameters occur, peaks are occasionally replaced by troughs 
in one parameter and not the other, and the amplitude of the 
oscillations in line centroid and width (and thus "total" inten- 
sity) oscillations are significantly reduced - these effects are 
also clearly present in the EIS data shown in Fig. [4] 

The results of the Monte Carlo simulations thus strongly 
indicate that the observed oscillations with EIS are compati- 
ble with a double component plasma where there is a quasi- 
periodically recurring upflowing component at 5% of the sta- 
tionary component's intensity, and with a relative velocity of 
order 50 km/s. This provides evidence for a scenario in which 
upflows cause some of the observed oscillations in intensity 
(e.g., in the SOHO, TRACE, STEREO, and EIS imaging data) 
and velocity and line width (in CDS and EIS spectra). Our 
results also indicate that the observed, less-than-ideal, cor- 
relations between the parameters are not a sign of a lack of 
correlation, but instead are exactly what is expected from the 
significant impact of finite spectral resolution, instrumental 
broadening and photon noise in the EIS instrument. We note 
that the line widths we have assumed (20-30 km/s) are sim- 
ilar to what is observed at the locations where EIS observes 
oscillations (after subtracting instrumental broadening). 

How unique is this interpretation? How well can we con- 
strain the properties of the second, faint, component? Are 
there other combinations of two Gaussian profiles (with one 
periodically recurring) that would lead to the same oscilla- 
tions in line intensity, centroid and width? To address this 
issue, we performed a parameter search in 4 dimensions for 
the case of profiles with S/N=30: varying a from values of 3 
to 9%, Sv from 10 to 80 km/s, and <tq and <j\ independently 
from 20 to 40 km/s. For each of the resultant timeseries of line 
(peak and total) intensity, centroid and width, we calculate the 
standard deviation over the course of one hour. We also cal- 
culate the ratio of the standard deviation of the peak intensity 
to that of the total intensity, and similar ratios for centroid to 
width and peak intensity to the centroid. We also calculate the 
standard deviation of the blueward asymmetry normalized to 
the peak intensity. For the sake of comparison we analyze the 
EIS observations of the locations at y = 60 — 65" and de- 
termine similar values for the observed oscillations. We use 
these observed standard deviations to constrain the parameter 
range of the forward models. This works surprisingly well, as 
shown in Fig. [7] By combining all 4 spectral line parameters 
with well-chosen ratios of those parameters, we find that for 
a given value of <tq and o\, only a small range of amplitudes 
a and velocity offsets Sv are compatible with the oscillations 
observed with EIS. Fig.[7]shows that for <tq = 30 km/s, and 
o\ = 20 km/s, only second components with amplitudes of 
order 4-8% and velocity offsets of order 30-50 km/s (blue- 
ward of the background component) are compatible with the 
observed behavior. This region is shown as the only region in 
the bottom right panel that is not hashed out. All other regions 
are excluded by constraints imposed by the EIS observations. 

As one might expect, these results depend significantly on 



the signal to noise assumed. For example, lower S/N observa- 
tions will naturally lead to relatively higher intensity pertur- 
bations as a direct consequence of photon noise. However, we 
have made sure that our assumed S/N of 30 matches that of the 
S/N of the observed constraints (derived from measurements 
based on summing over 2 spatial pixels, with each individual 
pixel having a S/N of ~20). 

The results also depend on the values assumed for ao and 
<7i. An accompanying movie (available in the online edition 
of the journal) illustrates that while the shape of the region 
of allowed parameters changes somewhat for different values 
of a, the range remains locked around amplitudes of order 
3-8% and velocities of 30-50 km/s. In fact, we find that for 
widths less than 20 km/s and more than 35 km/s, there are no 
solutions that are compatible with the observed oscillations 
(under our assumptions). Moreover, for realistic values of the 
line widths (which are of order 20-30 km/s after subtracting 
instrumental broadening, see also § 5), there are no solutions 
for offset velocities of order 10 km/s (which would be nec- 



essary for the wave interpretation of Wang et al. (2009b I, see 
§ 6). 

The forward model presented here assumes that the upflow- 
ing component is perfectly periodic, with identical amplitude 
throughout the timeseries. In reality, the upflowing compo- 
nent is likely only quasi-periodic, with an amplitude that can 
vary erratically. As a result, the standard deviation of each 
parameter might be different, which could impact the range 
of allowable parameters. However, the most stringent con- 
straints are actually the ratios of the standard deviation of var- 
ious parameters. These ratios are unlikely to change much if 
the upflowing component were variable in time and period, 
since all parameters will be impacted equitably by a lack of 
oscillatory power. Nevertheless, because of some of these 
limitations it is prudent not overinterpret the forward model. 
Consequently, we conclude that a second component that has 
a range (as opposed to exact values) of parameters (3-8%), 
and offset velocities (30-50 km/s) can explain the observed 
oscillations. 

In summary, we find that the results of the forward model 
strongly support the scenario posed in this Paper, that of 
a faint, but quasi-periodically recurring upflow component 
which impacts the intensity, centroid and line width of a sin- 
gle Gaussian fit to the spectral line profile. 

5. R-B "GUIDED" DOUBLE FITS 

Can we find direct evidence of this second component by 
performing double Gaussian fits? Motivated by these Monte 
Carlo simulations, and the presence of the line profile asym- 
metries in multiple spectral lines we now look to see if we 
can directly characterize the additional emission component. 
To do this we explicitly choose a double Gaussian component 
fit as one that is minimally consistent with the data (of course 
more free parameters will ensure a "better" fit to the data). In 
this case we chose not to perform a "blanket" double Gaussian 
profile fit, that is one where the fit parameters are completely 
free, or one where a physical premise is used to deduce, or 
loosely im pose, wher e the components lie in the spectral do- 
main (e.g., Peter 2000). We take advantage of the fact that the 
line profile asymmetries and large line widths are co-spatial 
with locations where the quality of fit measure, (i.e., the re- 
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duced x z ) is anomalously largaH We use the presence of a R- 

7 If the x 2 ma P that results from the single Gaussian fit to the line pro- 
files is not spatially "flat" then the spectra are not adequately described by a 
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Figure 8. Sample GA double Gaussian fits for the Hinode/ElS Spectroheli- 
ogram of Fig.Mjin the Fe XIV 274 A emission line. From top to bottom, left to 
right, the figure shows the peak intensity, Doppler velocity, and line width of 
the core component, the relative brightness, Doppler velocity and line width 
of the second component. 

B asymmetry in the line profile of a given pixel to determine 
that a double fit is needed, and to provide the fitting algorithm 
with starting parameters for that second component. For this 
first test of such a procedure we use a Genetic Algorithm (GA) 
Gaussian fitter t hat is based on a we ll tested, robust, although 
slow method ( |McIntosh et al.|1998) . 

This guided double fit procedure is quite simple. We find 
all locations where the R-B asymmetry exceeds a given mag- 
nitude (in this case 1% of the peak line intensity). We 
then approximate the R-B profile (as a function of veloc- 
ity/wavelength) by a Gaussian, and use the Gaussian only to 
estimate the center (in wavelength) of the R-B enhancement. 
The centroid of the initial line profile fit (used to compute the 
original R-B asymmetry) and the centroid of the R-B asym- 
metry profile which, depending on the sign of the R-B asym- 
metry will be on the blue (negative) or red (positive) side of 
the initial (core) profile, are the only values supplied to the GA 
Gaussian double fit. This algorithm then undertakes a global 
minimization of the double Gaussian line spectra fit allow- 
ing both core and second component centroids to move by 1 
spectral pixel to the blue or red of the supplied position. For 
reference, at 195 A in EIS this would allow a shift of ±35km/s 
from the input parameters for both emission components. We 
note also that all other parameters in the fit are free to span the 
range of possible values. 

Figure W shows the results of applying this technique to the 
spectro-heliogram shown in Fig.fl for the Fe XIV 274A emis- 
sion line. While rough, panels A, B, and C are consistent with 
their contemporaries in Fig. [T] while panels D, E and, F show 
the relative brightness (Jwing/Icore) of me second (wing) 
emission component, its velocity relative to the first compo- 
nent, and the line width respectively. The regions of large 
profile asymmetry under the line EIS slit for the timeseries 
observation show a second emission component of approxi- 
mate amplitude 6-8%, a velocity of 45-65km/s and (Gaussian) 
line width of 25- 35km/|D The widths of the second and first 
component are similar. 

single Gaussian; i.e., there is something m issing in the underlying physical 
description of the data (e.g., Bevington & Robinson 2003 1. 

8 The values for the width in this section exclude the instrumental broad- 
ening, which has been subtracted to allow comparison with the input values 
for the forward model. 
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Figure 9. Variation in the second GA double Gaussian fit components for 
the HinodefElS time series observation of Figs.|2|and plin the Fe XIV 264A 
emission line. The top panel shows the relative ^Brightness (Iwing/Icore) 
of the second (wing) emission component while the bottom panel shows its 
Doppler velocity (relative to the core component). For reference the dashed 
horizontal lines show the location of the time series chosen for Fig. [4] 

Fe XIV 264A: Comparison o1 R-B and Second Component 
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Figure 10. Extracted variation of the second Fe XIV 264A emission compo- 
nents of relative brightness (blue line) and velocity (green triangles) versus 
the variation in the (inverted) R-B diagnostic (red line, with blueward asym- 
metries positive). 

Figure|9]shows the results of applying the same R-B guided 
double Gaussian fit technique to the timeseries observations 
in the Fe XIV 264A emission line. This line is formed at 
the same temperature (in equilibrium) as Fe XIV 274A and 
is used here to illustrate that similar spectral asymmetries re- 
sulting from the second emission components occur across 
a range of temperatures. The two panels show the relative 
brightness (Iwing/ Icore> t0 P) °f me second (wing) compo- 
nent and its position relative to the rest position of the line 
(bottom). We see that the intensity of the second compo- 
nent modulates, but that the velocity of the second compo- 
nent stays relatively stable at 35-60 km/s. The (Gaussian) line 
width of the second component (not shown) is also stable in 
magnitude at 25-35 km/s. 

To illustrate these properties we show for Fe XIV 264A 
(Fig. 10 1 the variation in the relative brightness (blue line) 
and velocity (green triangles) of the second component, ex- 
tracted from the region inside the dashed lines (cf. Fig. [2] 
and|4]l, with the percentage change in R-B assymmetry at the 
same location. As for Fig. [4] the R-B timeseries is inverted 
such that excursions in the blue wing have positive values. 
We see that the blue wing enhancements in Fe XIV 264A are 
accompanied by a strong second emission component that has 
a fairly steady velocity at ^50km/s from the first (core) com- 
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ponent. This is consistent with what is observed for Fe XIV 
274A. We note also that the second component properties are 
remarkably consistent with those determined from the sim- 
ple forward model presented in the previous section. In fu- 
ture work we hope to expand and explore this encouraging 
"guided" double component fit technique to more data sets of 
interest, but that is beyond the scope of the present effort. 

6. DISCUSSION 

Before we can close our analysis there are a couple of out- 
standing questions. First, could the oscillations in line width 
(of order 1 .5 km/s) be caused by changes in the thermal broad- 
ening? In short, this is highly unlikely and for the following 
reason. In the locations we study in the current work, the 
overall line broadening (<7i/ e , the 1/e width) is of order 60 
km/s, whereas the thermal width of th ese lines is of order 
20 km/s. We know that the a±/ e — \J a 2 NT + af nst + o\ h , 
with <Ji ns t = 22mA the instrumental broadening, and oth = 
yj2kT/m with k the Boltzmann constant and m the ion mass. 
If we assume that the slight changes in line width Aa 1 u are 
caused only by thermal width changes, we can calculate the 
following: 



and: 



°th = °l/e - C 7 NT ~ °mst = 2fcT/m, (2) 



AT = ma 1 / e A<j 1/e /k (3) 



With the observed values of Aa\/ e ~ 1.5 km/s, we esti- 
mate that AT w 0.66 MK. Such a large change in temperature 
would significantly change the ionization balance, and essen- 
tially shift the temperature so far away from the peak of the 
contribution function for each ion, that very drastic changes in 
intensity would be expected. Because these are not observed, 
it is unlikely that the line width variations are predominantly 
caused by temperature changes associated with sound waves. 
This is further bolstered by the fact that the observed quasi- 
periodic variations are in phase for ions formed ac ross a wide 
range of temperatures (see e.g., Pascoe et al. 2007 for simula- 
tions of how temperature changes from sound waves can lead 
to strong phase shifts between oscillations of lines from ions 
in neighboring ionization stages). 

It seems from this analysis that the line widths changes 
cannot be caused by thermal changes associated with sound 
waves that propagate on the structure responsible for the dom- 
inant emission. Is it possible that sound waves propagating 
on a steadily blueshifted component (compared to the dom- 
inant, stationary, one) can explain the observations? Such a 
scenario could in principle cause linewidth changes associ- 
ated with blueward centroid shifts and intensity peaks. How- 
ever, it would require a secondary component that is persis- 
tently blueshifted compared to the dominant component for 
the whole 3 hours of the timeseries. We believe that our anal- 
ysis using both the R-B asymmetry variations and the double 
fits strongly implies the presence of a second component of 
emission that occurs at high velocity (i.e., it is associated with 
high velocity upflows). Our results thus pose a challenge to 
wave models to reproduce the observed amplitude of changes 
in line width, the clear R-B asymmetry variations, and the 
correlation of line widths with R-B signals at high velocities 
(>70km/s). Such models require sophisticated synthetic ob- 
servables and numerical models of the solar atmosphere (e.g., 
Pascoe et al.||2007] l and are beyond the scope of the current 
paper. 




Figure 11. Results of a Monte Carlo simulation illustrate that oscilla- 
tion finding software cannot distinguish quasi-periodic signals arising from 
stochastically driven upflows from those of propagating waves. Panel (a) 
shows a TRACE 171 A image where we added a series of simulated upflow 
events that originate at the footpoint of the cool loop emanating from the red 
cross. These upflow events are assigned the properties we expect for spicules 
heated to coronal temperatures: lifetimes of 100s, intensity of order 10% of 
the background loop structure, apparent propagation speed of 70 km/s (cho- 
sen to be along the negative y-axis in this case), and slight widening and 
fading. Panel (b) shows the timing of the randomly chosen occurrence of 
the upflow events at the location of the red cross in panel (a) (red symbols 
along the bottom of the plot), with the original timeseries at the location of 
the blue cross in panel (a) (blue dashed line), the simulated timeseries (full 
black line), and its running time difference (red dashed). Note the slight time 
delay (because of the finite propagation speed) between the spicule occur- 
rence at the red location, and the peaks in the timeseries at the blue location. 
Panel (c) shows the phase speed (in Mm/s) for regions where the wavetracker 
found significant propagating signal in the 3.5 mHz (250 s) passband with 
white contours for intensity of the TRACE 171 A images. Panel (d) shows the 
wavelet power for the timeseries of panel (b), i.e., for the blue cross location. 



The presence of a high-speed upflowing component is di- 
rectly shown by the results of our study of double Gaussian 
fits to the time series in the previous section. Generally dou- 
ble fits for line profiles with such a weak secondary com- 
ponent are poorly constrained and the results of the fit are 
strongly affected by the initial guesses for the propert ies (in- 
tensity, velo city, width) of the second component ( |McIntosh| 
|et al.|fT998] ). The result is often highly noisy and so addi- 
tional, strong physical or observational constraints must be 
applied - the approach attempted above directly uses observa- 
tions, and is a first step in that direction. It is pertinent here 
to note the perils of single Gaussian fits to data such as those 
discussed above - those where a second emitting component 
is visible, either in the spectral or imaging data. Our analysis 
shows that the results of a single Gaussian fit to such data can 
significantly impact the physical interpretation in a quite un- 
intentional way. It is possible that double fit analysis, driven 
from a knowledge of multiple emitting components in a res- 
olution element, like that demonstrated here, can help reduce 
that ambiguity and potential misinterpretation. 

The results discussed in the sections above provide a com- 
pelling scenario in which upflows at high speed (of order 
50-150 km/s ^] and low brightness (of order a few percent of 

9 The viewing angle will determine the amplitude of the line-of-sight ve- 
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the background brightness) are responsible for quasi-periodic 
changes in line intensity, velocity, width, and blueward asym- 
metries that propagate away from magnetic footpoints along 
coronal loops. We do not claim that that this scenario can ex- 
plain all reported observations of sound waves in the corona. 
After all, we are only observing one active region, and there 
are too many uncertainties (see below) at this point in time. 
However, our results strongly suggest that before the wave in- 
terpretation is applied to intensity oscillations, a careful anal- 
ysis of higher order moments of spectra is required. 

The presence of a second component that flows upward at 
speeds that are of order 50 km/s higher than the dominant 
plasma component begs the question how these two compo- 
nents are connected. Using the current data it is impossible to 
say whether the two components occur close to each other on 
neighboring field lines that are indistinguishable because of 
the spatio-temporal resolution of EIS which is modest com- 
pared to the fine-scale structuring and fast dynamics of the 
low solar atmosphere, or whether the enormous superposition 
of coronal lines along the line-of-sight samples two widely 
separated populations of plasma. 

In the following we briefly speculate on whether the "flows" 
scenario could explain some of the observed properties per- 
ceived to belong to sl ow magnetoacoustic waves (for a re- 
view, see [Nakariakov 2008)? As one example, we suggest 
that the flows scenario can naturally explain the lack of clean 
oscillations in plage-related loops, and the predominance 
of non-periodic propagating disturbances th at occur all over 
quiet Su n and plage-related loops (Mcintosh & De Pontieu 
2009a b ). We note that these non-periodic propagating distur- 
bances are identical in appearance (intensity, apparent propa- 
gation, lifetime) to t hose that are quasi-periodic and a bound 
in the literature (e.g.,[Schrijver et al.|1999| [Winebarge r et al. 
20021 |Sakao et al.||20071 |McIntosh et al.||2010a| |McIntosh 
& De Pontieu 2009b|ab. If the flows were driven stochasti- 



cally by an as yet unknown mechanism that affects the mag- 
netic field and occurs on dominant photospheric timescales 
(from p-mode oscillations at 4-10 minutes, and granulation 
on timescales of 5-15 minutes), one would not expect to see 
oscillations everywhere. In addition, the oscillations that are 
observed would not be expected to be very clean or signifi- 
cant. That is exactly what we see with TRACE and EIS. In 
addition, the timescale of quasi-periodicity would often vary 
strongly from location to location (observed with TRACE, see 



de Moortel et al. 2002b a ). In fact, detailed analysis of spectra 
in quiet Sun network regions shows that the blueward asym- 
metries do indeed sometimes recur quasi-periodically, but of- 
ten do not, which is fully compatible w ith the TRACE ob- 
servations (Mcintosh & De Pontieu|2009b| ). We find similar 
behavior in the EIS data analyzed here. Fig. [5] shows the re- 
sults of a wavelet analysis that shows quasi-periodicities of 
order 5-12 minutes for all derived properties of the line pro- 
files throughout the timeseries that can be re adily compared 
to Fig. 14 of |McIntosh & De Pontieu"! ( |2009a[ l. 

What drives these upflows? Give the proposed association 
with upflow events in the lower atmosphere (e.g., spicules; 
|McIntosh & De Pontieu||2009b| l, it is natural to assume a 
mechanism in which reconnection caused by photospheric 
motions which drive/alter the small-scale magnetic field in 



locity, determined from our RB, double fit and Monte Carlo analysis, and of 
the velocity in the plane of the sky, determined from space-time plots. The 
viewing angle likely changes from the footpoints to the upper parts of the 
loop. 



a stochastic manner. We should stress that the energy and 
mass release may not always lie in the chromosphere and 
appear as a spicule, but could well arise in the transition re- 
gion. Given the range of observed photospheric motions and 
the prevalent timescales, quasi-periodicity would occur natu- 
rally. This quasi-periodicity can be accentuated because the 
observed coronal "oscillations " are typically identified using 
a wavelet analysis ( e.g.. |Nakariakov & King 2007 Sych & 



Nakariakov 



2008 



To illustrate the ill-posedness of interpreting quasi- 
periodicities in the corona using imaging data alone, we have 
performed further forward modeling using Monte Carlo sim- 
ulations. We start with a TRACE 171 A image sequence and 
add an artificial signal that consists of a series of propagat- 
ing events with properties derived from the analysis presented 
above: low intensity (10% of background intensity), short 
lifetime at the footpoints (100s), high apparent propagation 
velocity in the plane of the sky (70km/s), and gradual fading 
and slight expansion in width (over e-folding distance of, re- 
spectively 15 and 20 Mm), as a result of decreasing density in 
the upflow events with height, e.g., because of fast expansion 
while the upflows merge into the pre-existing corona. 

To study the impact of a stochastic driver, we allow the 
"heating" events to occur randomly in time at the loop foot- 
point with a random uniform distribution given by the IDL 
function randomu.pro. On average, we introduce n events 
during the timeseries of duration AT, with AT/n of order 
the granular timescale (e.g., 300s). We then use two anal- 
ysis techniques that are commonly used to detect "oscilla- 
tions" in coronal imaging data: a wavelet analysis of the time- 
series of running differences (in this case with 120 s time dif- 
ference) in individual pixels (Torre nce & Compo|1998| ), and 
our "wavetracker" software which focuses on detec ting prop- 
agating signals in a set of narrow frequency ranges (Mcintosh 



|et al. ||2Q08p . An example of the analysis is given in Fig. 11 
which shows the original TRACE 171 A image, the simu- 
lated and original timeseries in one location, and the results of 
both the wavetracker and wavelet analysis. We can see that in 
this particular case the random bunching in time of individual 
propagating heating events leads to a clearly detectable quasi- 
periodicity in the wavelet analysis, and a nicely propagating 
signal (with 70 km/s) in the wavetracker analysis. However, 
no waves are present in the original data, and no waves have 
been introduced in the simulated data. Only heating events 
with the properties that we have deduced from the EIS data 
are included. This underscores the severe limitations of os- 
cillation finding techniques that are based on only intensity 
information, such as coronal images. 

To investigate how severe this issue is, we performed a large 
number of similar simulations with a range of characteristic 
timescales AT/n between 200 and 700s. This preliminary 
study shows that the combination of wavetracker and wavelet 
analysis of the resulting timeseries detects quasi-periodic sig- 
nals ("oscillations") in about 15-50% of the realizations, de- 
pending on the amplitude of the heating events, the signif- 
icance levels used in the wavelet analysis, and whether the 
running difference or original data is used for the wavelet 
analysis. This implies that the observed quasi-periodicities 

10 In the specific case of the EIS data, it is possible that the pointing jitter 
perpendicular to the slit (which cannot be corrected for and can be of order 
1-2") sometimes acts as a secondary cause of some of the quasi-periodicities 
in Figsp][T]and|2] We do not believe it is the dominant cause given the fact 
that sucn effects would impact the whole slit simi larly throughout th e whole 
timeseries. This is not what we observe (see also, |Wang et al.|2009b| 
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in the upflows are not necessarily a sign of a wave driver. 
More importantly, it illustrates that we cannot use imaging 
data alone to distinguish between upflow and wave scenarios, 
casting signficant doubt on the coronal seismology applica- 
tions of "slow-mode magneto-acoustic" waves. Our results 
suggest that such studies will require careful analysis of at 
least the first three moments of simultaneously observed spec- 
tral line profiles. 

Can the "flows" scenario explain what causes the rapid fad- 
ing and uniquely upward apparent motions of the propagating 
disturbances? In the wave scenari o, this has previously been 
attributed to damping of waves (Klimchuk 2006). In our sce- 
nario the visibility of the upward flows would critically de- 
pend on the density of the disturbances and on how the preex- 
isting coronal plasma absorbs these mass injections. Without 
a proper model of what drives these upflows, the associated 
heating and the impact on the corona, it is very difficult to 
predict how the density (and velocity!) of these upflow events 
would change with height. Perhaps the coronal part of the up- 
flow events tracks the behavior of the chromospheric counter- 
part (type II spicules) which show a clear decreas e of density 
with height in Hinode/SOT Ca II H timeseries (de Pontieu 
et al.||2007] )? In other words, perhaps the rapid decrease of 



ing disturbances of small amplitude observed on ubiquitous 
featur es observed in polar coronal holes, polar plumes (see, 
e-g 



Deforest & Gurman][l998) |Ofman et aij[l999l [Baner 



density with height in the upflow events is responsible for the 
apparent fading of the disturbances with distance along the 
coronal loops? Such a fading could perhaps also explain the 
fact that Wang et al. (2009b) found longer "periods" higher 
in the loop fan than at the footpoint. If only a fraction of the 
upflow events were dense enough to be visible at distances of 
20 Mm from the loop footpoints, the recurrence of the mo- 
tions observed higher in the corona would then tend to be 
on longer timescales, since some of the events recurring on 
shorter timescales would be invisible far away from the loop 
footpoints. 

An interesting twist to this scenario is added by preliminary 
analysis of AIA and Hinode data of upflow events, which in- 
dicates that the upflows are associated with heating of some 
rapidly moving chromospheric plasma to coronal tempera- 
tures. It is clear that such rapid heating could, in principle, 
lead to the rapid propagation into the corona of both ther- 
mal conduction fronts and/or sound waves. Depending on the 
initial conditions and details of the heating events, thermal 
conduction fro nts can pro pagate even faster than the speed 
of sound (e.g., [Hansteen 1993). Perhaps the strong heating 
and strong acceleration associated with these heating events 
then leads to the presence of both flows and associated ther- 
mal conduction front/sound waves? We believe that in the 
dataset considered here the flows dominate the observed os- 
cillations, but it is clear that careful analysis of spectra, and 
advanced numerical modeling will be required to disentangle 
which effects dominate in each observation. 

The impact of the results discussed in this Paper are wide- 
ranging. Coronal seismology based on the waves interpreta- 
tion of these propagating disturbances in coronal loops has 
been used to investigate the isothermal nature of loops, ther- 
mal substructuring of loops, the leakage of waves from the 
photosphere, the nature of thermal conductivity, etc (e.g., 



Robbr echt et al.||2001[ |de Moortel et al. 2002a b ; Kin g et al. 
2003 1 Ide Moortel & Rosner | |2007t |De Moortel & Bra dshaw 
20081 IManska et al.||2008| [Van Doorsselaere et aT|2 008b 
Ofman & Wang ||2008| |Owen et al l|2009| |De Moortel||2009, 
Wang et al. 2009a, De Moort er^Pascoe|2009[ |Marsh et at 
2009 1 |Marsh & Walsh 2009). Indeed, a similar slow-mode 



[jee et al.||2000| |2009| to name only a few of the many pa- 
pers on the topicX We believe that our results may affect 
the interpretation o f polar plume disturbances as well. |McIn-| 



tosh et al. (2010b I showed that the similarity in propagation 
speeds, amplitudes and quasi-periodicities of disturbances on 
polar plumes and those observed in other magnetic locations 
of the solar atmosphere suggests that the motions observed on 
plumes may not be the result of the passage of MHD waves, 
but a consequence of the quasi-periodic appearance of driven 
upflows. While we have not yet performed the kind of detailed 
spectroscopic analysis presented here on a polar plume (or an 
equivalent in an equatorial coronal hole, which would pro- 
vide the best opportunity for detecting profile asymmetries), 
it seems that the (almost global) properties and ubiquity of 
these upflow events indicates that a similar mechanism may 
be at play at the roots of the fast solar wind. This is illustrated 



MHD wave interpretation exists for high velocity propagat- 



by the SOHO coronal hole observation analysis of De Pontieu 
|et al.| ( [20Q9|) (and that of an upcoming paper - Mcintosh, Lea- 
mon & De Pontieu, "The Spectroscopic Footprint of the Fast 
Solar Wind", submitted to the ApJ). 

While a more extensive statistical study of oscillations us- 
ing spectroscopic data is required, it is clear that if our inter- 
pretation of the data presented is correct (i.e., that some of the 
disturbances are the result of driven upflows originating in the 
lower solar atmosphere), then much of the coronal seismology 
effort based on the interpretation of longitudinal, compres- 
sive, slow MHD modes in coronal regions not directly rooted 
in sunspots has to be reconsidered. The upflow interpretation 
on the other hand opens a new window into the connection 
between the chromosphere and corona with signi ficant impli- 
cation s for coronal heating mechanisms (e.g., |de Pontieu et al.| 
[2007] l. 

7. CONCLUSION 

We have demonstrated that recent EIS measurements of in- 
tensity and velocity oscillations, interpreted as direct evidence 
for propagating slow-mode waves in the corona, are actually 
accompanied by oscillations in the line width, and recurring 
asymmetries in line profiles across a range of temperatures. 
These facts, at the least, imply that coronal seismology using 
such propagating disturbances in the corona is not as straight- 
forward as assumed. Moreover, we show that these changes 
in intensity, velocity, line width and blueward asymmetry are 
compatible with a scenario in which faint upflows at high 
speed occur quasi-periodically and cause oscillations in the 
parameters determined from single Gaussian fits to spectral 
line profiles. Our results indicate that a significant fraction 
of the quasi-periodicities observed with coronal imagers and 
spectrographs that have previously been interpreted as prop- 
agating magnetoacoustic waves, may instead be caused by 
these quasi-periodic upflows. The uncertainty in the identi- 
fication of the physical cause for coronal oscillations signif- 
icantly impacts the prospects of successful coronal seismol- 
ogy using propagating, slow -mode magneto-acoustic waves. 
At the same time, the association of these propagating dis- 
turbances in coronal loops with high speed upflows provides 
an exciting new window into the mechanism that propels hot 
plasma into the corona. 
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Figure 1. Extracted EIS timeseries for the Fe XIII 202A (top), Fe XII195A (middle) and Fe XIV 274A (bottom) emission lines from 03:06-()6:00UT. This figure 
is similar to Fig.E] but for a different location (y=148"). In each panel we show the detrended percentage change in signal amplitude or velocity change in km/s. 
The green trace snows the temporal variance in the peak line intensity while the black, red, and blue traces show the changes in total profile intensity, Doppler 
velocity (line center shift), and line width respectively. Compare and contrast with the black dashed and dotted lines which show the 74-157km/s R-B analysis 
and detrended form (as for V and V nt ). The zero lines for these are drawn at -2 and -4 respectively and to stress our convention that values above that zero line 
represent an excess in the blue wing of the line. For clarity, the timeseries for intensity is offset by +4, velocity by +2 and the bottom linewidth by -2 from the 
zero-point on the y-axis. The presence of a second component in this location is also indicated in Fig. [9] 

APPENDIX 

In the appendix we show two more examples (Figs.[T]and[2]i of correlated oscillations in intensity, line centroid, linewidth, and 
blueward asymmetries. Here we did not sum over 5 pixels (as in Fig. |4]l, b ut over just two pixe ls. The locations where we have 
found the oscillations are in the same region as the two examples shown by Wang et al. (2009b ). 
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Figure 2. Extracted EIS timeseries for the Fe XIII 202A (top), Fe XII195A (middle) and Fe XIV 274A (bottom) emission lines from 03:06-06:00UT. This figure 
is similar to Fig.E] but for a different location (y=106"). In each panel we show the detrended percentage change in signal amplitude or velocity change in km/s. 
The green trace snows the temporal variance in the peak line intensity while the black, red, and blue traces show the changes in total profile intensity, Doppler 
velocity (line center shift), and line width respectively. Compare and contrast with the black dashed and dotted lines which show the 74-157km/s R-B analysis 
and detrended form (as for V and V nt ). The zero lines for these are drawn at -2 and -4 respectively and to stress our convention that values above that zero line 
represent an excess in the blue wing of the line. For clarity, the timeseries for intensity is offset by +4, velocity by +2 and the bottom linewidth by -2 from the 
zero-point on the y-axis. 



